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Abstract. A new weak Galcrkin (WG) finite element method is introduced and analyzed in this 
paper for the biharmonic equation in its primary form. This method is highly robust and flexible 
in the element construction by using discontinuous piecewise polynomials on general finite element 
partitions consisting of polygons or polyhedra of arbitrary shape. The resulting WG finite element 
formulation is symmetric, positive definite, and parameter-free. Optimal order error estimates in a 
discrete H 2 norm is established for the corresponding WG finite element solutions. Error estimates 
in the usual L 2 norm are also derived, yielding a sub-optimal order of convergence for the lowest 
order element and an optimal order of convergence for all high order of elements. Numerical results 
are presented to confirm the theory of convergence under suitable regularity assumptions. 
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1. Introduction. This paper is concerned with numerical methods for the bihar- 
monic equation with boundary conditions. For simplicity, consider the model problem 
that seeks an unknown function u — u{x) satisfying 

(1.1) A 2 u = f, into, 

(1.2) u = (, ondfl, 

d 

(1.3) ^r=4> onffi, 

on 

where A is the Laplacian operator, £ = C{ x ) an d 4> = 4>(x) are given functions defined 

on the boundary of the domain f2. Assume that the domain is open bounded with 
a Lipschitz continuous boundary <9f2 in R d , d = 2, 3. 

A natural variational formulation for the biharmonic equation (1.1) with Dirichlet 
and Neumann boundary conditions (1.2) and (1.3) seeks u € H 2 (Vl) satisfying u\gn — 
(, and an = 4> such that 

(1.4) (Au,Av) = (f,v), V«etf 2 (O), 

where Hq(£1) is a subspace of the Sobolev space 7J 2 (f2) consisting of functions with 
vanishing value and normal derivative on dfl. Based on the variational form (1.4), one 
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may design various conforming finite element schemes for (1.1)-(1.3) by constructing 
finite clement spaces as subspaces of H 2 (fl). It is known that _ff 2 -conforming methods 
essentially require C 1 -continuous piecewise polynomials on a prescribed finite clement 
partition, which imposes an enormous difficulty in practical computation. Due to 
the complexity in the construction of C 1 -continuous elements, ff 2 -conforming finite 
element methods are rarely used in practice for solving the biharmonic equation. 

As an alternative approach, nonconforming and discontinuous Galcrkin finite el- 
ement methods have been developed for solving the biharmonic equation over the 
last several decades. The Morley element [6] is a well-known example of nonconform- 
ing element for the biharmonic equation by using piecewise quadratic polynomials. 
Recently, a C° interior penalty method was studied in [2, 3]. In [8], a hp- version 
interior-penalty discontinuous Galerkin method was developed for the biharmonic 
equation. To avoid the use of C 1 -elements, mixed methods have been developed for 
the biharmonic equation by reducing the fourth order problem to a system of two 
second order equations [I, 4, 5, 7, 10]. 

In this paper, we will develop a highly flexible and robust weak Galcrkin finite 
clement method for the biharmonic equation that allows the use of generalized ap- 
proximating functions on general partitions consisting of polygons or polyhedra of 
arbitrary shape with certain shape regularity. The weak Galcrkin method refers to 
a numerical technique for partial differential equations where differential operators 
are interpreted and approximated as distributions over a set of generalized functions. 
The method/idea was first introduced in [11] for second order elliptic equations, and 
the concept was further developed in [12, 13, 9]. By design, weak Galerkin methods 
use generalized and/or discontinuous approximating functions on general meshes to 
overcome the barrier in the construction of "smooth" finite element functions. 

Intuitively, a weak Galerkin finite element scheme for the biharmonic equation 
(1.1)-(1.3) can be derived by replacing the differential operator A in (1.4) by a discrete 
weak Laplacian, denoted by A w . However, such a straight forward replacement may 
not work without including a mechanism that enforces a certain weak continuity of the 
underlying approximating functions. A weak continuity can be realized by introducing 
an appropriately defined stabilizer, denoted as s(-, •). Formally, our weak Galerkin 
finite element method for (1.1)-(1.3) can be described by seeking a finite element 
function Uh satisfying 

(1.5) {A w u h ,A w v) + s(u h ,v) = (f,v) 

for all testing functions v. The goal of the paper is to specify all the details for (1.5), 
and further justify the rigorousness of the method by establishing a mathematical 
convergence theory. 

The paper is organized as follows. In Section 2, we introduce some standard 
notations for Sobolev spaces. A weak Laplacian operator and its discrete version will 
be introduced in Section 3. In Section 4, we shall present two WG finite element 
schemes for the biharmonic equation (1.1)-(1.3). In Section 5, we shall introduce 
some local L 2 projection operators and then derive some approximation properties 
which arc useful in a convergence analysis. In Section 6, we shall establish optimal 
order error estimates for the WG finite element approximation in a H 2 -equivalent 
discrete norm. In Section 7, we shall derive an error estimate for the WG-FEM 
approximation in the usual L 2 norm. Results from two numerical experiments are 
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reported in Section 8. Finally, we provide some technical results in the appendix that 
are critical in dealing with finite element functions on arbitrary polygons/poly hedra. 

2. Preliminaries and Notations. Let D be any open bounded domain with 
Lipschitz continuous boundary in R d , d = 2,3. We use the standard definition for the 
Sobolcv space H S (D) and the associated inner product (-,-) s .d, norm || • || Si d, and 
seminorm | • \ s ,d for any s > 0. For example, for any integer s > 0, the scminorm 
| • \ aD is given by 



\v\., D = ( e jHsm 



K \a\=s " 



with the usual notation 



a = (ai,...,a d ), |a| = ai + . . . + a d , d Q = JJd";. 

3 = 1 



The Sobolev norm || • \\ m ,D is given by 



\V\\m,D 



The space H°(D) coincides with L 2 (D), for which the norm and the inner product 
are denoted by ||-||_d and (•, respectively. When D = 0, we shall drop the subscript 
in the norm and inner product notation. 

The space 7J(div; D) is defined as the set of vector-valued functions on D which, 
together with their divergence, are square integrable; i.e., 

H(div; D) = {v : v £ [L 2 (D)] d , V • v G L 2 (D)} . 

The norm in iJ(div; D) is defined by 




3. Weak Laplacian and Discrete Weak Laplacian. For the biharmonic 
problem (1.1)-(1.3) with the variational form (1.4), the principle differential oper- 
ator is the Laplacian A. Thus, we shall introduce a weak Laplacian operator defined 
on a class of discontinuous functions. For numerical purposes, we will define a discrete 
version of the weak Laplacian in polynomial subspaces. 

Let T be any polygonal or polyhedral domain with boundary dT. By a weak 
function on the region T we mean a function v = {v 07 Vb 7 v g } such that v £ L 2 (T), 
vi, £ Hi(dT), and v g ■ n £ H~i(dT), where n is the outward normal direction of T 
on its boundary. The first and second components vo and Vb can be understood as the 
value of v in the interior and on the boundary of T. The third component v 9 intends 
to represent the gradient of v on the boundary of T. Note that vt, and v g may not be 
necessarily related to the trace of vq and V«o on dT, respectively. 



4 



Denote by W(T) the space of all weak functions on T; i.e., 

(3.1) W(T) = {v = {v ,v b ,w g } : v G £ 2 (T), v b G #*(0T), v s • n G fl-*(0T)}. 
Let (-, -)qt be the inner product in L 2 (dT). Define G 2 (T) by 

G 2 (T) = {^: ^^(nA^i 2 ^)}. 
For any ip G G 2 (T), we have G H{div 1 T) 1 and hence • n G H~i(dT). 

Definition 3.1. TTie dwaZ of L 2 (T) can be identified with itself by using the 
standard L 2 inner product as the action of linear functionals. With a similar inter- 
pretation, for any v G W{T), the weak Laplacian of v = {v n ,v b} v g } is defined as a 
linear functional A w v in the dual space o/G 2 (T) whose action on each ip G G 2 (T) is 
given by 

(3.2) (A w v, ip) T = (v , Aip) T - (v b , Vip ■ n) aT + (v g • n, ip) dT , 

where n is the outward normal direction to dT. 

The Sobolev space H 2 (T) can be embedded into the space W(T) by an inclusion 
map i w ■■ H 2 (T) -> W(T) defined as follows 

iw{4>) = WW, <t>\ 9T , V0| aT }, G H 2 (T). 

With the help of the inclusion map iw, the Sobolev space H 2 (T) can be viewed as a 
subspace of W(T) by identifying each <fi G H 2 (T) with iw{4>)- Analogously, a weak 
function v — {v a , v b , v g } G W(T) is said to be in H 2 (T) if it can be identified with a 
function <p G H 2 (T) through the above inclusion map. It is not hard to see that the 
weak Laplacian is identical with the strong Laplacian in H 2 (T); i.e., A w v = Av for 
all functions v G H 2 (T). 

Next, we introduce a discrete weak Laplacian operator by approximating A w 
in a polynomial subspace of the dual of G 2 (T). To this end, for any non- negative 
integer r > 0, denote by P r {T) the set of polynomials on T with degree no more than 
r. A discrete weak Laplacian operator, denoted by A mjr j, is defined as the unique 
polynomial A WiTt TV G P r {T) satisfying the following equation 

(3.3) (A Wir , T v,ip) T = (v , Aip) T - (v b , V(p ■ n) 9T + (v 9 • n,cp) 9T , VV G P r {T). 

4. Weak Galerkin Finite Element Schemes. Let Th be a partition of the 
domain f2 into polygons in 2D or polyhedra in 3D. Assume that Th is shape regular in 
the sense as defined in [12] (seel also Appendix A). Denote by £h the set of all edges 
or fiat faces in Th, and let £\ — £h\dH be the set of all interior edges or fiat faces. 

For any given integer k > 2, denote by Wk(T) the discrete weak function space 
given by 

(4.1) W k (T) = {{v ,v b ,v g } : v G P k (T), v b G P k (e), v 9 G [P k -i{e)] d , e C dT}. 

By patching W k (T) over all the elements T G Th through a common value on the 
interface £®, we arrive at a weak finite element space Vh given by 

V h = {{vo,v b ,v g }: {vo,v b ,v g }\ T G W k (T), VT G %}■ 



■5 



Denote by A/> the trace of Vh on dil from the component v b . It is clear that 
consists of piecewise polynomials of degree k. Similarly, denote by T/j the trace of Vh 
from the normal component of v 9 as piecewise polynomials of degree k — 1. Denote 
by V h ° the subspace of Vh with vanishing traces; i.e., 

V£ = {v = {v ,v b ,v g } eV h ,v b \ e = 0, v 3 -n| e =0, ecdrndo}. 

Denote by A Wt k-2 the discrete weak Laplacian operator on the finite element 
space Vh computed by using (3.3) on each element T for k > 2; i.e., 

(A Wik -2v)\ T = A Wik -2,T{v\ T ), Vu G V h - 

For simplicity, we shall drop the subscript fc — 2 in the notation A Wy k-2 for the discrete 
weak Laplacian. We also introduce the following notation 

(A w v,A w w) h := {A w v,A w w) T - 

Ten 

4.1. Algorithm I. For any Uh — {u n ,u b ,u g } and v — {v ,v b ,v g } in Vh, we 
introduce a bilinear form as follows 

s(u h ,v) := ^2 h^CVuo - u g , Vv - w g ) dT + ^ h T 3 (u Q - u b , v Q - v h ) dT . 
Ten Ten 

Weak Galerkin Algorithm 1. A numerical approximation for (l.l)-(l.S) can 
be obtained by seeking Uh — {u , u b , u g } e Vh satisfying u b — Q b ( and u g ■ n = Q gn <fr 
on dQ and the following equation: 

(4.2) (A w u h ,A w v) h + s(u h ,v) = (f,v ), V v = {v ,v b ,v g } G V h °, 

where Q b ( is the standard L 2 projection onto the trace space Ah and Q gn <p is the L 2 
projection onto the normal component of the gradient trace space Th- 

The following is a useful observation concerning the finite element space V h °. 
Lemma 4.1. For any v e V h a , let |||i>||| be given as follows 

(4.3) IHII 2 = (A w v, A w v) h + s(v, v). 

Then, \\\ ■ \\\ defines a norm in the linear space V h ° . 

Proof. We shall only verify the positivity property for ||| • ||| . To this end, assume 
that I I = for some {vo 7 v b} v g } G V h °. It follows that A w v = 0, v = v b , and 
V«o = v s on each element T or dT, as appropriate. We claim that Av = holds 
true locally on each element T. To this end, for any ip G Pu-2{T) we use A w v = 
and the definition (3.3) to obtain 

(4.4) = (A w v,<p) T 

= (vo, A(p) T - (v b , V<p ■ n) dT + (v 9 • n, <p) dT 

= (Au , <p)t + (vo - v b , Vip ■ n) dT + (v g ■ n - Vu • n, ip) dT 

(4.5) = (Avq, <p)t, 
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where we have used the fact that v — v b = and Vv — v g = in the last equality. 
The identity (4.5) implies that At> — holds true locally on each element T. This, 
together with v = Vb and V«o = v g on 9T, shows that v is a smooth harmonic 
function globally on f2. The boundary condition of Vb = then implies that heOoii 
ft, which completes the proof. □ 

Lemma 4.2. The weak Galerkin finite element scheme (4-2) has a unique solu- 
tion. 

Proof. Let and be two solutions of the weak Galerkin finite element 
scheme (4.2). It is clear that the difference e/j = — u^ 2 ' is a finite element 
function in V h ° satisfying 

(4.6) (A w e h ,A w v) h + s{e h ,v) =0, V v = {v , v b , v 9 } e Vg. 
By setting v = e h in (4.6) we obtain 

(A w e h , A w e h ) h + s(e h ,e h ) = 0. 
It follows from Lemma 4.1 that e/j = 0. This shows that = u h 2 \ □ 

4.2. Algorithm II. Here we describe another weak Galerkin finite element 
scheme that has less number of unknowns than (4.2). This second WG-FEM scheme is 
formulated in a subspace of Vh, denoted by Vh, that uses only the normal component 
of the vector v s for any v — {v , Vb, v g } G Vh- To be more precise, let us introduce a 
set of normal directions on Eh as follows 

(4.7) T>h = {n e : n e is unit and normal to e, e G £h}- 
The weak Galerkin finite element space Vh is given as follows 

(4.8) V h = {v = {v 0l v b ,v g n e } : v G Pk(T),v b G Pk(e),v g G Pn(e),eC dT}, 

where v g can be viewed as an approximation of Vv ■ n e . Denote by V® the subspace 
of Vh with vanishing traces; i.e., 

V° = {v = {v ,v h ,v g n e } G V h ,v b \ e = 0, v g \ e = 0, e C dTndfl}. 

For any Uh = {u$, Ub, u g n e } and v = {vq, Vb, v g ii e } in Vh, we introduce a bilinear 
form as follows 

s(u h ,v) := ^ h^CVuo ■ n e - u g , Vv ■ n e - v g ) 9T + ^ h T 3 (u - u b , v - v b )dT- 
Ten Ten 

Weak Galerkin Algorithm 2. A numerical approximation for (l.l)-(l.S) can 
be obtained by seeking Uh = {uq, u b , u g n e } G Vh satisfying Ub = Q b ( and u g = Q gn <t> 
on d£l and the following equation: 

(4.9) (A w Uh,A w v) h + s(uh,v) = (f,v Q ), V v = {v ,Vb,v g ii e } e V%, 

where Qb( is the standard L 2 projection onto the trace space Ah and Q gn <f> is the L 2 
projection onto the normal component of the gradient trace space Th- 

Like (4.2), the WG-FEM scheme (4.9) can be proved to have one and only one 
solution in the corresponding finite element space. Details are left to interested readers 
for a verification. 
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5. L 2 Projections and Approximation Properties. For each element T, 
denote by Qo the L 2 projection onto P^{T),k > 2. For each edge/face e C dT, 
denote by Q b and Q s the L 2 projection onto Pfc(e) and [Pt_i(e)] d , respectively. Now 
for any u e H 2 (il), we can define a projection into the finite element space Vu such 
that on the element T 

QhU = {Qqu, Q b u, Qg(Vu)}. 

In addition, denote by Qh the local L 2 projection onto Pk~2(T). 

Lemma 5.1. On each element T G Th, the L 2 projections Qh and Qh satisfy the 
following commutative property with the Laplacian A and the discrete weak Laplacian 

(5.1) A w (Q h u) = Q h (Au) 

for all u £ H 2 (T). 

Proof. For any t e Pk-2(T), it is not hard to see that 

(A w Q h u, t)t = {Qqu, At) t + (Q g (Vu) • n, r) dT - (Q b u, Vr • n) dT 
= (u, At) t + (Vu • n, t) 9t - (u, Vt • n) dT 
= (Au,t) t 
= (QhAu, t )t, 

which implies the desired identity (5.1). □ 

The commutative property (5.1) indicates that the discrete weak Laplacian of the 
L 2 projection of smooth functions is a good approximation of the Laplacian of the 
function itself in the classical sense. This is a nice and useful property of the discrete 
weak Laplacian in application to algorithm design and analysis. 

The following lemma provides some approximation properties for the projection 
operators Qh and Qh- 

Lemma 5.2. Let Th be a finite element partition of Q satisfying the shape regu- 
larity assumption as specified in [12]. Then, for any < s < 2 and 2 < m < k we 
have 

(5.2) Y, h 2 T s \\u-Q u\\l T <Ch 2 ^\\u\\ 2 m+1 , 
Ten 

(5.3) ]T h 2 T s \\Au-® h Au\\l T < Ch 2 ^\\u\\ 2 n+1 . 

Ten 

Proof. The proof of this lemma is similar to that of Lemma 5.1 in [12], and the 
details are thus omitted. □ 

Using the trace inequality (A. 3) in the Appendix with p = 2, we can derive the 
following estimates which are useful in the convergence analysis for the weak Galerkin 
finite element schemes (4.2) and (4.9). 



8 



Lemma 5.3. Let 2 < m < k, w e i2" max { m+1 ' 4 }(Q) ; and v e V h . There exists a 
constant C such that the following estimates hold true: 



(5.4) MAw-QfcAHlir) <Ch m - 1 \\w\\ m+1 , 
\TeT h ) 

(5.5) ( Y 4llv(Aw; - q h A W )\\ 2 dT j < c/^aiHIm+i + foWHU). 

\TeTh / 

(5.6) ( ^ h^\\V(Q w) - Q fl (Vtu)Hlr j < CTi^lHlm+i, 
\TeTh / 

(5.7) ( ^ ^HQo^-QHUt) ^C^'IML+i- 

Here Sij is the usual Kronecker's delta with value 1 when i = j and value otherwise. 

Proof. To derive (5.4), we use the trace inequality (A. 3) and the estimate (5.3) 
to obtain 



Y h T \\Aw - Q h Awf dT 



Ten 



<C J2 {\\Aw-® h Aw\\ 2 T + h 2 T \\V(A w -® h Aw)\\ 2 T ) 



Ten 

n-2\\ || 

m+1- 



< Ch 2m - 2 \\w" 2 



As to (5.5), we use the trace inequality (A. 3) and the estimate (5.3) to obtain 

Y h 3 T \\V(Aw-Q h Aw)\\ 2 T 
Ten 

<C Y (h 2 T \\V(Aw - Q h Aw)\\ 2 T + h 4 T \\^ 2 (Aw - Q h Aw)\\ 2 T ) 
Ten 

<Ch 2m - 2 (\\w\\ 2 m+1 + h 2 S m<2 \H\ 2 ). 

As to (5.6), we have from the definition of Q g , the trace inequality (A. 3), and the 
estimate (5.2) that 

Y h T 1 \\V(Q w)-Q g (Vw)\\ 2 T 
Ten 

< Y h^WVQow- Vw\\ 2 T 
Ten 

<C Y (h T 2 \\VQ w - Vw\\ 2 T + \\VQ w - VHIi,t) 
Ten 

<Ch 2m - 2 \\w\\ 2 m+1 . 
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Finally, we use the definition of Qb and the trace inequality (A. 3) to obtain 
^2 h T 3 \\Qo w ~ Qbw\\ 2 dT < ^2 h T 3 \\Qo w - w\\It 

TeT h TeT h 

<C (h T 4 \\Qow-w\\ 2 T + h T 2 \\W(Q w-w)\\ 2 T ) 
Ten 

<Ch 2m - 2 \\w\\ 2 m+1 . 
This completes the proof of (5.7), and hence the lemma. □ 

6. An Error Estimate in H 2 . The goal here is to establish an error estimate 
for the WG-FEM solution Uh arising from (4.2) and (4.9). For simplicity, we will 
focus on the error analysis for (4.2) only; the analysis can be easily modified to cover 
the WG-FEM scheme (4.9) without any difficulty. 

First of all, let us derive an error equation for the weak Galcrkin finite element 
solution. 

Lemma 6.1. Let u and Uh = {u ,Ub,u g } e Vh be the solution of (l.l)-(l.S) and 
(4-2), respectively. Denote by 

eh = QhU - u h 

the error function between the L 2 projection of u and its weak Galerkin finite element 
approximation. Then the error function e/j satisfies the following equation 

(A w e h , A w v) h + s(e h ,v) = ^ (Au - Q h Au, (Vv - v s ) • n) dT 

TeT h 

(6.1) - ( V ( Au - ^hAu) -n,v - v b ) dT + s(Q h u, v) 

Ten 

for all v e V° . 

Proof. Using (3.3), the integration by parts, and the fact that A w QhU = Qh(Au) 
we obtain 

(A w Q h u, A w v) T 

= (v , A(A w Q h u)) T + (v 9 • n, A w Q h u) dT - (v b , W{A w Q h u) ■ n) 9T 
= (Av , A w Q h u) T + (v , \7(A w Q h u) ■ n) dT - (Vv ■ n, A w Q h u) dT 

+ (v g • n, A w Q h u) dT - (v b , V(A w Q h u) ■ n) dT 
= {Av , A w Q h u) T + (v - v b , V{A w Q h u) ■ n) dT - ((Vv - v s ) • n, A w Q h u) dT 
= (Av ,Q h Au) T + (v - v b , V(Q h Au) ■ n) dT - {(Vv - v s ) • n, Q h Au) dT 
= (Au, Av ) T + (v - vb, V(Q h Au) ■ n) dT - ((Vv - v g ) ■ n, Q h Au) 9T , 

which implies that 

(Au, Av ) T = (A w Q h u, A w v) T - (v - v b , V(Q h Au) ■ n) dT 

(6.2) + {(Vv -v g )-n,Q h Au) dT . 

Next, it follows from the integration by parts that 

(Au, Av )t = (A 2 u, v ) T + (Au, Vv ■ n) 9T - (V(Au) • n, v ) dT . 
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Summing over all T and then using the identity (A 2 u, v ) = (/, v a ) we arrive at 

53 ( A «, Aw )t = (/, vo) + 53 ^ Au ' ( Vw o ~ v s) ' n ^T ~ (V(Att) ■ n, u - «b)dT, 
Ter h TeTft Ten 

where we have used the fact that v s • n and v b vanishes on the boundary of the domain. 
Combining the above equation with (6.2) leads to 

(A w Q h u, A w v) h = (/, v ) + 5Z (Am - Q h Au, (Vw - v s ) • n) aT 

TeT h 

- 53 (V(Au - Q h Au) -n,v - v b ) dT . 
Ten 

Adding s(QhU, v) to both sides of the above equation gives 

{A w Q h u, A w v) h + s(Q h u, v) = (/, v ) + 53 ( Au _ QhAu, (Vv - v g ) ■ n) 9T 

Ten 

(6.3) - ^ (V(Au - Q h Au) ■ n, v - v b ) dT + s(Q h u, v). 

Ten 

Subtracting (4.2) from (6.3) yields the following error equation 

{A w e h , A w v) h + s(e h ,v) = 53 ( Au - QhAu, (Vw - v fl ) • n) dT 

Ten 



- 5Z _ QhAu) -n,v - v b ) aT + s(Q h u, v) 



Ten 



h 



for all v G Vfi- This completes the derivation of (6.1). □ 

The following is an estimate for the error function eh in the trip-bar norm which 
is essentially an iJ 2 -norm in Vfi. 

Theorem 6.2. Let Uh G Vh be the weak Galerkin finite element solution arising 
from (4-2) with finite element functions of order k > 2. Assume that the exact solution 
of (l.l)-(l.S ) is sufficiently regular such that u G iJ max { fe+1 ' 4 }(fi). Then, there exists 
a constant C such that 

(6.4) \\\u h - Q h u\\\ < Gh k ~ x (|M| fe+ i + W fe ,2||u|| 4 ) • 

In other words, we have the optimal order of convergence in the H 2 norm. 

Proof. By letting v = eh in the error equation (6.1), we obtain the following 
identity 

(6.5) \\\e h \l = J2(Au- Q h Au, (Ve - e s ) • n) aT 

Ten 

- 53 ( V ( Au - QhAu) ■ n, e - e b ) dT 
Ten 

+ 5^ ^ T 1 ( V( 9ow - Q 9 Vu, Ve - e g ) dT 



Ten 

+ 51 h T 3 (Q u - Q b u, e - e b ) dT . 
Ten 
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Now using the Cauchy-Schwarz inequality and the estimates (5.4) and (5.5) of Lemma 
5.3 with m = k one arrives at 



(6.6) 



( Au - QhAu, (Ve - e fl ) • n) 9T 



Ten 



< { E M Au - QfcA«|||r ) ( E ^'HVeo - e g \\ 2 dT 

\TET h 



\TeT h 



|u||fc+i|||eft|| 



and 
(6.7) 



• n, e - e b )dT 



Ten 



< E 4I|V(A U - QfcAuJUl,, X ^ 3 lko - e 6 | 
\Ter h / \Ter h 

^C^-HllulU+i + ^Mllull^lefcll. 



Similarly, it follows from the Cauchy-Schwarz and the estimates (5.6) and (5.7) that 



(6.8) 

and 
(6.9) 



X hj, 1 (VQqu - QgVu, Ve -e g )dT 
Ten 



< Ch 



k-l 



u\\k+i\\\eh\\ 



X h T 3 (Q u-Q b u, e ~e b )dT 
Ten 



< Ch 



fc-ii 



|w||fe+i|||eh|||. 



Substituting (6.6)-(6.9) into (6.5) yields 

||e h ||| 2 <Cft fc - 1 (|H| fc+ i + M fc , 2 ||«||4) III HI , 
which implies (6.4). This completes the proof of the theorem. □ 

7. An Error Estimate in L? . This section shall establish an estimate for the 
first component of the error function in the standard L 2 norm. To this end, we 
consider the following dual problem 



(7.1) 
(7.2) 
(7.3) 



A 2 w = eo in fi, 
w = 0, on dSl, 
Vw ■ n = on dil. 



Assume that the above dual problem has the following regularity estimate 
(7.4) H|4<C||eo||. 

Theorem 7.1. Let Uh € Vh be the weak Galerkin finite element solution arising 
from (4-2) with finite element functions of order k > 2. Let to = min(fc,3). Assume 
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that the exact solution of (l.l)-(l.S ) is sufficiently regular such that u G H k+1 (Q) 
and the dual problem (7.1)-(7.3) has the H 4 regularity. Then, there exists a constant 
C such that 

(7.5) HQou - «o|| < ^ fe+t °- 2 (|| U || fe+1 + M M || U || 4 ). 

In other words, for quadratic elements (i.e., k = 2) we have a sub-optimal order of 
convergence given by \\Qou — uo\\ = 0(h 2 ). But for cubic or higher order of elements, 
we have the optimal order of convergence \\Qou — tto|| — 0(h k+1 ). 

Proof. Testing (7.1) with the error function eo on each element and then using 
integration by parts to obtain 

|| eo || 2 = (AVeo) 

= ^2(Aw,Ae ) T + Y (V(Aio) • n,e ) aT - Y (Aw, Ve • n) aT 

TeT h TeT h TeT h 

= Y ( Aw ' Ae o)r + Y (V(Aio) ■ n, e - e b ) dT - Y ( Aw > ( Ve o - e s) ' n )dT> 
Ter h Ter h Ter h 

where we have used the fact that e b and e s • n vanishes on the boundary of il. Using 
(6.2) with w in the place of u, we can rewrite the above identity as follows 

||e || 2 = (A w Q h w, A w e h ) h + Y (( v ( Aw ) - V{Q h Aw)) ■ n, e - e b ) dT 

TeTh 

- Y ( Aw ~ ^Aw, (Ve - e g ) ■ n) aT . 

TeT h 

Next, we have from the error equation (6.1) that 

{A w Q h w, A w e h )h = Y ( Au ~ ^ hAu > ( V( 3o^ - Q 9 Vw) • n) 9T 

TeT h 

- Y ( v ( Aw ~ QhAu) ■ n, Q w - Q b w) 9T 
Ten 

- s(e h , Q h w) + s(Q h u, Q h w). 
Combining the above two equations we obtain 

(7.6) ||e || 2 - Yl <( V ( Aw ) - V(QhAw)) ■ n, e - e b ) dT 

TeT h 

- Y ~~ ® h AW ' ( Ve ° ~~ ^ ' U ^dT ~ S ( e h, Qhw) 

Ten 

+ Y ( Au - ^Au, {VQow - Q g Vw) ■ n) dT 

Ten 

- Y ( v ( Am - QhAu) ■ n, Q w - Q b w) dT + s(Q h u, Q h w). 
TeTh 

Each of the six terms on the right-hand side of (7.6) can be bounded by using the 
Cauchy-Schwarz inequality and Lemma 5.3 as follows. 
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that 
(7.7) 



For the first term, it follows from the estimate (5.5) and the fact t = min(3, k) < 3 
<V(Au> - Q h Aw) ■ n, e - e b )dT\ 



Ten 



<( h 3 T \WAw-® h Aw)\\ 2 dT ) ( J2 h T 3 \\eo-e b \\ 2 dT 



\Ten 

<Cft*°- 1 (||t«|| t0+ i + /wJ t0 ,2||«;||4)||e fc | 
HUIIM!- 



\TeT h 



<C/i*°- 11 



As to the second term, we use the Cauchy-Schwarz inequality and the estimate 
(5.4) with m = t to obtain 



(7.8) 



(Aw - QhAw, (Ve - e g ) ■ n) dT 



Ten 



< h T \\Aw-Q h Aw\\ 2 aT ) ]T ^HVeo-eJir 



\Ten 
<Ch t °- 1 \ 



\Ten 



w \\to+l\\\ e h\\ 



KCh^WwMehl 
Analogously, for the third term, we have 
(7.9) \8(e h ,Q h w)\<Ch t °- 1 \\w\\4e h l 



The fourth term can be bounded by using Lemma 5.2 as follows 



(7.10) 



( Au - QhAu, (VQ w - Q 9 Vw) ■ n) aT 



Ten 



< Y, h\\^-Q h Au\\ 2 aT \ J2 h-'WVQow-VwWlA 
\TeT h ) \Ten ) 

<c^- 1 || u || fe+1 ^- 1 IMk +1 
<c^+ to - 2 || u || fe+1 || w || 4 . 

As to the fifth term, we again use the Cauchy-Schwarz inequality and Lemma 5.3 
to obtain 



(7.11) 



Y (V(A« - Q h Au) • n, Q w - Q b w) dT 
Ten 



< E h s T \\W(Au-q h Au)\\ 2 dT \ ( ]T h T 3 \\Q w-w\\ 2 aT 
\Ten / 

< Ch k -\\\u\\ k+1 + hS^WuW^h^WwW 

< Ch k+t °- 2 (\\u\\ k+1 + h6 k , 2 \\u\\ 4 )\\w\\ 4 . 



\Ten 

to + 1 
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The last term on the right-hand side of (7.6) can be estimated as follows. 
(7.12) \s{Q h u,Q h w)\ 



< 



h^iVQou - Q g \7u,\7Q w - Q g Vw) aT 



Ten 



^ h T 3 {Q u - Q b u, Q w - Q b w) dT 



Ten 



< h^WVQou- Vu| 



dT 



\Ten 



\TeT h 



dT 



+ h T 3 \\Qou-u\\ 2 dT j j 

< C^-^luHfc+i ^-'H^llto+i 
<Ch k+t °- 2 \\u\\ k+1 \\w\U. 



fr T 3 iiQo?« - ^1 



dT 



Substituting the estimates (7.7)-(7.12) into (7.6) yields 

II e 1| 2 <C^- 1 ||e h |||||t«||4 + C'/i fc+t0 - 2 (||«|| fc+ i+Wfc,2ll«ll4)l|t«||4. 

Using the regularity estimate (7.4) we arrive at 

||e || < Ch^WlehlW + Ch k+t "- 2 (\\u\\ k+1 + hS k , 2 \\u\U). 

Finally, by combining the above estimate with the H 2 error estimate (6.4) we obtain 
the desired L? error estimate (7.5). □ 

8. Numerical Results. Our numerical experiments are conducted for the weak 
Galerkin finite element scheme (4.9) by using the following finite element space 

V h = {v = {v ,v b ,v g ii e }, v G P 2 (T), v b 6 P 2 (e), v g e Pi(e),T e T hl e G £ h }. 

For any given v = {vq, v bl v g n e } e Vh, the discrete weak Laplacian, A w v, is computed 
locally on each element T as a function in Pq(T) by solving the following equation 

(A w v, VOt = (vo, A^)t + (v g n e ■ n, V)st - (v b , Vip ■ n) 9T , 

for all ip e Po{T). Since V S -Po(r), the above equation can be simplified as 

(8.1) {A w v, ip) T = (v g n e ■ n, ^) aT . 

The error for the WG-FEM solution is measured in two norms defined as follows: 



(8.2) IK 1 2 :=J2 ( I ' \A w v h \ 2 dx + h^ f \(\7v )-n e -v 
+/iy 3 / (v - v b ) 2 ds ) , 



l ds 



dT 



(8.3) |K|| 2 := V f \v \ 2 dx, 

rn^T- JT 



(A discrete ff 2 -norm), 
(Element-based L 2 -norm). 



Ten- 
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8.1. Test Case 1. Here we consider the forth order problem that seeks an un- 
known function u = u(x, y) satisfying 

A 2 u = f 

in the square domain — (0, l) 2 with homogeneous Dirichlet and Neumann boundary 
conditions. The exact solution is given by u = x 2 (l — x) 2 y 2 (l — y) 2 , and the function 
/ = f(x,y) is computed to match the exact solution. Uniform triangular meshes 
are used and they are constructed as follows: (1) partition the domain into n x n 
sub-rectangles; (2) divide each square element into two triangles by the diagonal line 
with a negative slope. The mesh size is denoted by h = l/n. Table 8.1 shows that 
the convergence rates for the WG-FEM solution in the H 2 and L 2 norms are of order 
0(h) and 0(h 2 ), respectively. 



Table 8.1 

Test Case 1: Numerical error and convergence rates in H 2 and L 2 . 



h 


\\u h - Qhu\\ 


order 


||«o - Qou\\ 


order 


2.5000e-01 


2.5683c-01 




3.3304c-02 




1.2500e-01 


1.3540c-01 


9.2359e-01 


9.1046c-03 


1.8710 


6.2500c-02 


7.2378c-02 


9.0360c-01 


2.6049c-03 


1.8054 


3.1250c-02 


3.8275e-02 


9.1915e-01 


7.3257c-04 


1.8302 


1.5625e-02 


1.9687e-02 


9.5916e-01 


1.9461e-04 


1.9124 


7.8125c-03 


9.9457c-03 


9.8510e-01 


4.9762c-05 


1.9675 



8.2. Test Case 2. Here we solve the biharmonic equation on the domain of unit 
square with nonhomogeneous Dirichlet and Neumann boundary conditions. The exact 
solution is given by u — sin(7ra;) sin(7rj/), and the function / = ,f(x,y) is computed 
accordingly. 

The mesh triangulations are constructed in the same way as in the test case 1. 
The numerical results are presented in Table 8.2 which confirm the theory developed 
in earlier sections. 

Table 8.2 

Test Case 2: Numerical error and convergence rates for the biharmonic equation with non- 
homogeneous boundary conditions. 



h 


\\u h - Qhu\\ 


order 


||m - Qou\\ 


order 


2.5000c-01 


2.4536e+01 




3.1862c+00 




1.2500e-01 


1.2794c+01 


9.3943e-01 


8.5298c-01 


1.9013 


6.2500e-02 


6.7243c+00 


9.2801c-01 


2.3439c-01 


1.8636 


3.1250c-02 


3.4811e+00 


9.4984e-01 


6.2578c-02 


1.9052 


1.5625c-02 


1.7657e+00 


9.7930e-01 


1.6066c-02 


1.9616 


7.8125e-03 


8.8709e-01 


9.9309e-01 


4.0534c-03 


1.9868 



More numerical experiments should be conducted for the WG-FEM schemes (4.2) 
and (4.9). In particular, it would be interesting to see results for WG approximations 
of high order elements. It is also important to hybridize (4.2) and (4.9) so that 
unknowns related to vq can be eliminated locally on each element. 
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Appendix A. The goal of this Appendix is to establish some fundamental esti- 
mates useful in the error estimate for general weak Galcrkin finite element methods. 
First, we derive a trace inequality for functions defined on the finite element parti- 
tion Th with properties as specified in [12]. For completeness, we review the shape 
regularity assumption in the following subsection. 

A.l. Domain Partition and Shape Regularity. Let Th be a partition of the 
domain f2 consisting of polygons in two dimensions or polyhedra in three dimensions 
satisfying a set of conditions to be specified. Denote by Eh the set of all edges or flat 
faces in Th, and let £° = Eh\dVL be the set of all interior edges or flat faces. For every 
element T € Th, we denote by |T| the area or volume of T and by hr its diameter. 
Similarly, we denote by |e| the length or area of e and by h e the diameter of edge or 
flat face e e Eh- We also set as usual the mesh size of Th by 

h = max fix- 

Ten 




Fig. A.l. Depiction of a shape-regular polygonal element ABCDEFA. 

All the elements of Th are assumed to be closed and simply connected polygons or 
polyhedra, see Fig. A.l. The partition Th is said to be shape regular if the following 
properties are satisfied. 

Al: Assume that there exist two positive constants g v and g e such that for every 
element T e Th we have 

(A.l) Qv h d T < Qeht 1 < \e\ 

for all edges or flat faces of T. 
A2: Assume that there exists a positive constant k such that for every element 
T G Th we have 

(A.2) Kh T < h e 

for all edges or flat faces e of T. 
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A3: Assume that the mesh edges or faces are fiat. We further assume that for every 
T £ Th, and for every edge/face e £ dT, there exists a pyramid P{e,T, A e ) 
contained in T such that its base is identical with e, its apex is A e £ T, and 
its height is proportional to hr with a proportionality constant a e bounded 
from below by a fixed positive number a*. In other words, the height of 
the pyramid is given by a Jit such that er e > a* > 0. The pyramid is also 
assumed to stand up above the base e in the sense that the angle between 
the vector x e — A e , for any x e £ e, and the outward normal direction of e is 
strictly acute by falling into an interval [0,9o] with # < §• 

A4: Assume that each T £ Th has a circumscribed simplex S(T) that is shape 
regular and has a diameter hsrr) proportional to the diameter of T; i.e., 
hs(T) < l*hr with a constant 7, independent of T. Furthermore, assume 
that each circumscribed simplex S(T) intersects with only a fixed and small 
number of such simplices for all other elements T £Th- 

A. 2. A Trace Inequality. The following is trace inequality concerning func- 
tions in the Sobolev space W 1,p . 

Lemma A.l (Trace Inequality). Let Th be a partition of the domain VL into poly- 
gons in 2D or polyhedra in 3D. Assume that the partition Th satisfies the assumptions 
(Al), (A2), and (A3) as specified in Section [12]. Let p > 1 be any real number. 
Then, there exists a constant C such that for any T £ Th and edge/face e £ dT , we 
have 

(A.3) \\e\\i P{e) < chS (\\e\\ 

Lp(T) 

where 9 £ W 1 ^^) is any function. 

Proof. We shall provide a proof for the case of e being a flat face and £ C 1 (T). 
To this end, let the flat face e be represented by the following parametric equation: 

(A.4) x e = 77) :=(M£,v),M£,v),'h&v)) 

for (£, rf) £ D C M 2 . By Assumption A3, there exists a pyramid P(e, T, A e ) with apex 
A e = x* := (xl,X2,xl) contained in the element T. This pyramid has the following 
parametric representation: 

(A.5) x(t,£,»7) = (l-i)0(£,»7)+tx, 

for (i, £, i]) £ [0, 1] x D. For any given x e £ e, the line segment joining x e and the 
apex can be represented by 

x(i) = x e + t(x* - x e ). 
From the fundamental theorem of Calculus, we have 

\0\r (x e ) - |0|P(x(i)) = - f d T (\6\r(x e + TL0))dT, u = x, - x e . 
Jo 

The above can be further rewritten as 

|0|f(xe) - |C(x(t)) = -p f W-^g^O) (W • to) dr. 

Jo 
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Let q = -^j. It follows from the Cauchy-Schwarz inequality that for any t G [0, |] we 
have 



|0|*(xe) < |0|*(x(t)) +p (j 2 W\w\dr\ " (j 2 |V0|* 



CO 



Using the Young's inequality ab < edP + e pb q for any non-negative real numbers 
a, b, and e > 0, we arrive at 

i_ i 

\6\ p (x e )<\8\ p (x(t))+ P e-i j 2 \6\ p \u;\dT + p e [ 2 \V9\ p \lo\(1t. 

Jo Jo 

Now we integrate the above inequality over the flat face e by using the parametric 
equation (A. 4), yielding 

f |0| p (x e )|<^ x ^dr, < [ \9\ p (x(t))\<t>i x ^dr, 

JD J D 



(A.6) +pe—p \6\ p \uj\\^ x ^Id^dr 

Jo Jd 

(A.7) + P e f 2 f \W6\ p \oj\\^ x <f>„\dSdridT. 

Jo Jd 

Observe that the integral of a function over the following prismatoid 
Pi :={x(t,£,»7): (t, £, rf) € [0, 1/2] x D} 

is given by 

f /(x)dx= / 2 [ f{ X ( T ,t,ri))J(T,£,ri)d£dridT, 

JPi JO JD 

2 

where J(r, £, 77) = (1 — r) 2 |(^ x 0,,) • w| is the Jacobian from the coordinate change. 
The vector <j>£ x <f) v is normal to the face e, and oj = x* — x e is a vector from the base 
point x e to the apex x*. The angle assumption (see Assumption A3 of Section ??) 
for the prism P(e, T, A e ) indicates that the Jacobian satisfies the following relation 

(A.8) J(n?,i|)>j|^xA,| H, re [0,1/2] 

for some fixed /j,q G (0,1). Observe that the left-hand side of (A.6) is the surface 
integral of \8\ p over the face e. Thus, substituting (A.8) into (A.6) and (A.7) yields 



\8\ p de< j^\6\ p (x(t))\<j> ( x^dr, 

+4^ V? f |6fdx + 4 w - 1 e f \V9\ p dx. 
Jp L Jpi 

2 2 

Now we integrate the above with respect to t in the interval [0, \] to obtain 

(A.9) ^Jj8\ p de<JjJ D \e\ p ( X (t))\^x^d V dt 

+2p/i - 1 £-" f \9\ p dx + 2pfi 1 e [ \V9\ p dx. 
JPi Jpi 
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Again, by substituting (A. 8) into the right-hand side of (A. 9) we arrive at 

(A.10) lfj 6 \ Pde ^ ^o 1 ^ 1 " 1 / 2 j D \V\ P (x(t))J(t,!i,v)d!id V dt 

+2p^ l e~i [ \e\ p dx + 2 P fi 1 e [ \V9\ p dx. 
JPi JPi 

The first integral on the right-hand side of (A. 10) is the integral of \9\ p on the pris- 
matoid Pi . It can be seen from the Assumption A3 that 

(A.ll) < a*h T 



for some positive constant a*. By taking e — l? T 1 , we have e p = hj} . It then 
follows from (A.10) and (A.ll) that 

\8\ p de < Ch T x ^ J ^ \6\ p dx + h r J ^ |V6»| p dxj , 

which completes the proof of the Lemma. □ 

A. 3. A Domain Inverse Inequality. Next, we would like to establish an esti- 
mate for the L p norm of polynomial functions by their L p norm on a subdomain. To 
this end, we first derive a result of similar nature for general functions in W 1,p . 

Lemma A. 2. Let K C M. d be convex and v E W^ P {K) with p > 1. Then, 
(A.12) M\l P(K) < 2 -^\\v\\ P LP(s) + {^P-J WW 



IIP 

ILp(K)' 



where 5 is the diameter of K , S is any measurable subset of K , and lu^ — ^(l/l) * s 
the volume of unit ball in R d . 

Proof. Since C 1 (if) is dense in W 1,P (K), it is sufficient to establish (A.12) for 
v e C 1 (iir). For any x,y £ K, we have 

\v\ p (x) = \v\ p (y) - /' X ^ d r (\v\ p (x + ruj))dr, co = f^-. 

Jo |y~ x l 

From the usual chain rule and the Cauchy-Schwarz inequality we obtain 

•l x -y| 

(A.13) < \v\ p (y) +p e -f / y \v\ p dr+pe f V \d r v\ p dr, 



r\*--y\ 

M P (x) — \v\ p (y) — p / \v\ p ~ 1 sgn(v)d r v(x + ruj)dr 
Jo 



where e > is any constant. Let 

V(x) = 



o Jo 



\v\P{x), x6if 
0, x g K 
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r i^i^x), xeif 

W \ 0, xgA". 
Then, the inequality (A. 13) can be rewritten as 

_i r°° r 00 

M P ( X ) < M P (y) + p £ p / y(x + ruj)dr + p e I W(x + no)dr. 

Jo Jo 

Integrating the above inequality with respect to y in S yields 

(A.14) |S|M P (x) < / \v\ p dS 
Js 

POO poo \ 

J V(x + ruj)dr + e J W(x + rui)dr) dy. 



+p [e p 

J\x-y\<8 V 

It is not hard to see that 

p poo poo p pS 

/ / V(x + ruj)drdy = / / / V '(x + to) p d ~ x dp duj dr 

J\x-y\<sJo Jo J\ui\ = lJo 

fid poo P 

= — / / V(x + ruj)duj dr 
d Jo J\w\=i 

(A.15) =^|x-y| 1 - d H P (y)dy. 
Analogously, we have 

(A.16) / / W(x + rio)drdy = 6 - f |x - y| 1 - d |S r u(y)| p dy. 

J\x-y\<5 Jo a Jk 

Substituting (A.15) and (A.16) into (A.14) yields 

\s\ i«r(x) 

\v\ p dS + ^- (e~p |x - yr- d M P (y)dy + £ |x - yr- d |V«(y)| p dy 
Observe that the following holds true 

/ Ix-yl^dx^aSd-i, 

where S^-i is the "area" of the unit surface in R d . The volume of the unit ball in R d , 
LUd, is related to Sd-i as follows 

Sd-i 
d 

Now integrating both sides with respect to x in K gives 

\S\ [ \v\ p dK<\K\ f \v\ p dS + P LJ d S d+1 ( e -f f \v\ p dK + e f \Vv\ p dK) , 
Jk Js V Jk Jk J 

which yields the desired estimate (A. 12) by setting e — (^^Js\~^ ■ ^ 



Consider a case of Lemma A. 2 in which the convex domain K is a shape regular 
d-simplex. Denote by hx the diameter of K. The shape regularity implies that 



21 



1. the measure of K is proportional to h K , 

2. there exists an inscribed ball B K C K with diameter proportional to hx ■ 

Now let 5* be a ball inside of K with radius r$ > s*Iik- Then, there exists a fixed 
constant k* such that 

(A.17) \K\ < k.\S\. 

Apply (A.17) in (A. 12) and notice that t*\S\ > h d K and S = h K . Thus, 
(A.18) \H P LV[K) < 2k4v\\ p lp{S) + (2pnh K u d ) p HVwIl^^j. 

For simplicity of notation, we shall rewrite (A.18) in the following form 
(A.19) \\vf KtP < a \\vr Sp + ai hF K \\Vvr Ktp . 

If v is infinitely smooth, then a recursive use of the estimate (A.19) yields the following 
result 

n 

(A.20) \\v\f Kp <J2 a J h ^ jv \\l P + a n + ihT^ n+ ^ P K, P - 

j=o 

In particular, if v is a polynomial of degree n, then 

n 

(A.21) ll«IISr, P <E°^ll VJ ' v lls, P - 

3=0 

The standard inverse inequality implies that 

\\V j v\\s,p<hK i \\v\\ S j,. 
Substituting the above into (A.21) gives 

(A.22) \H P K , p <\\v\\l p . 
The result is summarized as follows. 

Lemma A. 3 (Domain Inverse Inequality). Let K C K d be a d-simplex which has 
diameter hx and is shape regular. Assume that S is a ball in K with diameter r$ 
proportional to hx; i.e., r$ > <;*hx with a fixed > 0. Then, there exists a constant 
C = C(<r*,n) such that 

(A.23) \H LP{ K)<C(^,n)\\v\\ LP(s) 
for any polynomial v of degree no more than n. 

A. 4. Inverse Inequalities. The usual inverse inequality in finite element anal- 
ysis also holds true for piecewise polynomials defined on the finite element partition 
Th provided that it satisfies the assumptions A1-A4. 

Lemma A. 4. Let Th be a finite element partition of ft consisting of polygons or 
polyhedra. Assume that Th satisfies all the assumptions A1-A4 and p > 1 be any real 
number. Then, there exists a constant C — C(n) such that 

(A.24) || V^|| T ,p < C(n)^ 1 |M| T , p , VT e % 
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for any piecewise polynomial ip of degree n onTh- 

Proof. The proof is merely a combination of Lemma A. 3 and the standard inverse 
inequality on d-simplices. To this end, for any T G 7ft, let S(T) be the circumscribed 
simplex that is shape regular. It follows from the standard inverse inequality that 

l P\\s(T),p- 

Then we use the estimate (A. 23), with K — S(T), to obtain 

||v^||t, p < ch^WvWsv < ch^y\\ T<P , 

where S is a ball inside of T with a diameter proportional to Iit- This completes the 
proof of the lemma. □ 

Lemma A. 5. Let Th be a finite element partition of £1 consisting of polygons or 
polyhedra. Assume that Th satisfies all the assumptions A1-A4 and p > r > 1 be any 
two real numbers. Then, there exists a constant C = C(n) such that 

(A.25) IMUf(n) < Chp~l \\f\\ L -(n) 

for any piecewise polynomial ip of degree n onTh- 

Proof. For any T G Th, let S(T) be the circumscribed simplex that is shape 
regular. It follows from the standard inverse inequality that 

d _ d 

IMkp < |MU(T),p < C/l£ r |M|s(T),r~ 

We then use the estimate (A. 23), with K — S(T), to obtain 

d_d 

II v\\ t, p < Ch£ r |M|s> 
(A.26) < C7if~'IMkr, 

where S is a ball inside of T with a diameter proportional to hr- It follows from 
(A.26) that 



MIl»(ii) = Mr* 



Ten 

< Ch*-* £ |M|?. r . 
Ten 

Since *j < 1, then we have from the above inequality that 



MIL>(0) < Chf- d MT,r 



\Ten 



<Chf- d Y,M r T,r, 

Ten 

which implies the desired inverse inequality (A.25). □ 
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